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Abstract 


_j)lhis  report  describes  a  general  purpose  computer  package  for  analysis 
of  the  composite  criterion  model  of  Srinivasan  and  Shocker#  (Puycliumeferika^ 
VOl.  3&t~  No .  A  |  pt  Ari) .  A  standard  format  allows  data  for  a  wide  range 
of  applications  to  be  analyzed  using  this  package.  The  report  describes 
the  development  of  the  package#  including  the  data  base  and  the  FORTRAN 
codes  for  performing  the  analysis.  The  sample  problem  used  by  Srinivasan 
and  Shocker  is  analyzed  to  demonstrate  the  package. 
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A  COMPUTER  PACKAGE  FOR  THE  COMPOSITE  CRITERION  MODEL 


The  composite  criterion  model  was  proposed  by  Srinivasan  and  Shocker 
[A]  as  a  technique  for  analyzing  a  decision  maker's  preferences  when  there 
are  multiple  conflicting  criteria.  The  model  is  based  on  the  assumed 
existence  of  a  value  function  which  assigns  a  cardinal  score  to  any  alter¬ 
native.  Associated  with  alternative  i  is  a  vector  [a^ ]  of  attribute 
values.  The  value  function  is  assumed  to  be: 


•i  '  vaay>)  -  |  V« 


(i) 


The  goal  of  the  composite  criterion  model  is  to  use  a  limited  number  of 
explicit  pairwise  comparisons  to  reveal  the  proper  values  of  the  weights, 
Wj,  for  a  given  decision  maker. 

Given  two  alternatives,  i  and  j,  if  the  decision  maker  prefers  i  over 
J,  then  it  is  assumed  that  s^  >  s^,  i.e., 

I - 1  Vjk : 0 

or  £  (*irV  Kk  £ 0  (2) 


Thus,  if  we  have  a  set,  fi,  of  pairwise  comparisons,  there  is  a  corresponding 
set  of  constraints,  (2),  which  must  be  satisfied  by  the  weights  w^. 

Since  it  is  possible  that  the  decision  maker  may  err  in  the  preference 
responses  (or  that  the  linear  value  function  assumption  is  not  completely 
satisfied),  the  feasible  set  defined  by  (2)  for  (i,  j)  e  ft  may  be  empty. 

In  that  case,  we  would  like  to  determine  a  "best"  set  of  weights. 

The  composite  criterion  model,  as  proposed  by  Srinivasan  and  Shocker, 
determines  the  set  of  weights  which  minimizes  the  sum  of  the  Infeasibilities 
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in  the  constraints  (2) : 


minimize  T  . 

(l.j)efl  iJ 

(3) 

jp  aijkWk  +  Zij  -  0 

(i,j)en 

(4) 

kcP  (i, j)sft  ijk 

v  c  1 

k 

(5) 

z±j  >° 

(i*j)  e  S 

(6) 

w,  >  0  k 

k  - 

ePl 

(7) 

v,  <  0  k 

k  - 

6  P2 

(8) 

w,  urs  k 

k 

EP3 

(9) 

p  *  pi  u  p2  u  *3; 

p±  n  Pj  -  0  i  1 5 

Note  that  in  order  to  determine  the  weights,  wfc,  the  linear  programming 
problem  (3)  -  (9)  must  be  solved. 

Any  composite  criterion  problem  is  defined  by  the  sets  0,  P^,  Pj,  P^ 
and  the  vectors  [a^^].  Thus,  it  is  reasonable  to  develop  a  general  pur¬ 
pose  data  base  and  solution  procedure.  The  package  we  have  developed  is 
based  on  the  overall  design  illustrated  in  Figure  1.  The  LP  optimizer 
is  the  XMP  suite  of  FORTRAN  routines  developed  by  Marsten  [l] .  The  entire 
package  is  written  in  FORTRAN  (ANSI  66),  is  reasonably  portable,  and  can 
accomodate  problems  with  as  many  as  30  restricted  attributes  (or  15  unre¬ 
stricted  attributes)  and  100  pairwise  comparisons. 

Parameter  sizes  on  the  CCMODL  program  can  be  increased  to  accommodate 


any  size  problem  with  the  amount  of  array  space  increasing  as  a  function 
of  the  model  parameters.  Let: 


MAXFPC  *  maximum  number  of  forced-paired  comparisons, 

MAXWTS  *  maximum  number  of  weights  that  can  be  specified 
where  every  unrestricted  weight  counts  as  two 
weights, 

A  ■  number  of  words  required  for  array  storage. 

Then: 

A  >  63  MAXFPC  +  7  MAXWTS  +  4  (MAXFPC)  (MAXWTS) 

12  (MAXFPC  +  1)  (MAXWTS)  (MAXFPC)  +  151 
+  MAXWTS  +  3  MAXFPC  +  1 

STANDARD  CCM  DATA  BASS 


The  data  base  contains  a  complete  description  of  the  problem,  from 


which  the  interface  routines  generate  the  LP  coefficients.  There  are  five 

record  groups  in  the  data  base: 

group  1:  this  is  a  title  record  to  allow  a  description  to 
appear  on  the  report;  (10A8)  format, 

group  2:  this  gives  the  magnitude  of  the  problem; 

(KATTA,  KST1ML1 ,  KFCOMP) ,  where, 

KATTA  =  number  of  attributes 
KSTMLI  =  number  of  stimuli  or  alternatives 
KFCOMP  «=  number  of  paired  comparisons 
(313)  format 


group  3:  a  set  of  flags  to  indicate  sign  restrictions  on 
the  weights;  the  first  record  is 


KWFLAG 


(12)  format; 


’1  if  P  - 
1-1  if  P  *  P2 
2  if  P  -  P3 
•0  otherwise 


If  KWFLAG  -  0,  then  KATTA  flags  KW(i)  are  required, 
one  for  each  attribute 


KW(i) 
(12)  format. 


1 

2 

3 


s 


group  4:  the  set  0;  KFCOMP  pairs  (1,  j)  are  required 
where  i  is  preferred  to  J 
(214)  format 

group  5:  attribute  values;  for  each  stimulus,  a  set  of 
KATTA  values  are  required,  in  (8F10.0)  format; 

KSTMLI  sets  are  required. 

Example .  Srinivasan  and  Shocker  describe  a  small  problem  with  five 
stimuli  and  three  attributes.  The  data  for  this  problem  is  given  below: 

«  *»  { (1,2) ,  (3,1),  (1,4),  (1,5),  (2,3),  (2,4),  (5.2),  (3,4),  (3,5),  (5.4)} 


*  P  -  P. 


Stimulus  Attribute  Vector 

1  (2,6,6) 

2  (3,4,7) 

3  (5,4,4) 

4  (6,3,2) 

5  (8,1,3) 

The  standard  CCM  data  set  for  this  problem  is  given  in  Table  1. 

*  There  is  a  small  difference  between  this  problem  and  the  Srinivasan  and 
Shocker  problem  in  [4],  the  difference  being  that  in  [4]  all  weights 
are  restricted  to  be  non-negative  while  in  this  problem  all  weights 
are  unrestricted. 


LP  OPTIMIZER 

The  LP  optimizer  is  the  XMP  package  developed  by  Marsten  [1].  XMP 
Implements  a  revised  simplex  algorithm  and  uses  the  LA05  subroutines  to 
manage  an  LU  factorization  of  the  basis  matrix  from  the  Harwell  library. 
XMP  is  written  in  a  ’/Vanilla”  FORTRAN  where  only  very  simple  features  ..  . 
of  the  language  have  been  used.  XMP  operates  in  main  memory  where  the 
memory  space  that  XMP  uses  must  »e  provided  by  the  calling  program.  XMP 
is  hierarchically  structured  and  any  XMP  subroutine  may  be  called  at  any 
level  of  the  hierarchy. 


The  design  objecCives  of  XMP  were  flexibility,  extendability,  and 
reliability  with  execution  time  as  a  secondary  objective.  XMP  contains 
46  subroutines,  10500  lines  of  FORTRAN  and  requires  10395  words  for  the 
object  code. 

INTERFACE  ROUTINES 

The  interface  code  MCCMDL,  consists  of  551  lines  of  FORTRAN  and  is 
internally  documented.  The  code  requires  1246  words  for  the  object  code 
and  29708  words  for  arrays  (this  includes  the  arrays  required  for  XMP). 

As  presently  configured,  MCCMDL  can  accommodate  100  stimuli,  30  attributes 
and  100  paired  comparisons.  The  number  of  paired  comparisons  is  the  cri¬ 
tical  parameter,  since  it  increases  both  the  number  of  variables  and  the 
number  of  constraints  in  the  LP  problem. 

The  CCMODL  program  was  used  to  solve  the  problem  described  in  the 
data  set  listed  in  Table  1,  the  output  is  shown  in  Table  2.  The  results 
shown  in  Table  2  are  the  same  results  obtained  by  Srinivasan  and  Shocker 
in  [4].  The  run  was  completed  in  .486  CP  seconds. 

As  discussed  in  [4]  and  as  seen  in  Table  2,  stimuli  1,  2,  5  have  the 
same  score.  This  occurred  because  of  the  number  of  intransitivities  in 
the  pairwise  comparisons  in  the  example.  Since  every  Z ^  except  ”  .02174, 
is  equal  to  zero  this  implies  that  every  preferred  stimulus  in  a  pairwise 
comparison  - has  a  higher  score  than  the  other  stimuli  in  the  comparison 
except  in  comparison  5.  In  comparison  5  the  Decision  Maker  preferred  sti-r 
mulus  2  to  stimulus  3  when  comparing  scores  we  see  that  stimulus  3  has  a  - 
higher  score  (by  .02174)  than  does  stimulus  2. 


A  listing  of  the  FORTRAN  source  for  MCCMDL  is  given  in  Appendix  A 
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Table  2.  Result  for  Srinivasan-Shocker  Example 


TRIAL  i  FCK  THE  C  .C  .  NO  DEL  FECGRAN,  ALL  WEIGHTS  ARE  UNREST  RICTEQ 


THE  OPTIMAL  CSJECTI VE  FUNCTION  VALUE  IS  Z=  .0217  4 


WEIGHT# 


LUE** 


1 

2 

3 


.2  3913 
.19565 
.15217 


2 


***  VALUE'**’ 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0 .0GG00 
0 .00000 
G . C  CO  00 
0  .C00C0 
.02174 
C  .00000 
0 .00000 
0 .00000 
0.00000 
G.  0  000  0 


STINILUS 


1 

2 

3 

4 

5 


**-SCORE** 


2.56522 
2.56522 
2 . 56696 
2.32609 
2 .56522 


PhOUHAfi  r;CCMDL(  INPUT .OUTPUT , TAPt3= INPUT ,TAPEb=GOTPUT ,TAPb5 ) 


C 

C  Thb  PURPOSE  Of  This  PROGRAM  IS  70  SOLVE  ThE  COMPOSITE  CRITERION  MODEL 

C  bl  ThE  USE  Of  LINEAR  PROGRAMMING.  ThE  L.P.  MODEL  WILL  bE  SOLVED  bY  THE 

C  AMP  SUBROUTINE  PACKAGs,. 

C 

C  SlhCE  ThE  Xi.P  PACKAGE  Of  SUonGUTINtS  IS  SO  WELL  DOCUMENTED  ThE  COMMENTS 

C  II.  THIS  PROGRAM  WILL  ONLY  fcE  DIElCTlD  TOWARD  ThE  VAHIABLES  AND  LOGIC  Of 

C  THIS  PROGRAM  ONLY.  TnE  PARTS  OF  THIS  PROGRAM  COVERED  li«  ThE  XhP  DGCUKENT- 

C  ATION  WILL  Lb  SET  Of  BY  TnE  FOLLOWING - , 

C 

C  This  PROGRAM  IS  WRITTEN  hi  RCr.ALD  K.  RUNYAN  UNDER  THE  DIRECTION  OF 
C  Dh.  LEON  F.  MCGINNIS. 

C 

C  GEGOKIA  Ir.STITUTE  OF  TECHNOLOGY,  PALL  19bO. 

C 

C  DIMENSION  STATEMENTS 

c 

DIMENSION  STIKLI ( 1 CO , 30 ) , KW ( 30 ) , W EIUhT( 30 ) , RHALF( 1 0 1 , 30  )  ,  Z2  ( 1 00  ) 

DIMENSION  HEADER ( 1 C ) , KWDLF ( 30 ) 

INTEGER  OMEGA (100,2) 

C 

C  FOR  ChANGES  IN  ThE  PARAMETERS, THE  ARRAY  SIZES  ShOULD  Eb  ThE 
C  FOLLOWING  SIZE. 

C  STIMLKMAXf  PC  .KAXV.TS) 

C  Kli(KAXViTS) 

C  W EIGHT (KAXWTS) 

C  RhALF(MAXFPC+1 .MAXWTS) 

C  Z2(MAXFPC) 

C  HEADER(IO) 

C  OMEGA ( NAXFPC , 2 ) 

C  KN’DEf  (  MAXWTS  ) 

C 
C 

C  DIMENSION  AND  VARIABLE  DEFINITION  FOP.  XMP 

REAL  b( 101 ) ,BASCB( 101 ) ,BASLB( 101 ) ,BASUB( 101 ) , 

X  B£TAR(101 ) , BOUND, CAhDA( 101 ,5) ,CANDCJ(5) ,CJ,COLA( 101), 

X  LJ, MEMORY (20700 ) ,UJ,UZEnO( 101 ) ,XB2ER0( 101 ) ,YQ( 101 ) ,Z 
C 

INTEGER  BNDTYP  ,COLL£N , COLMAX ,  DFEASQ ,  DTEhM , DUNBR , 

X  FACTOR, IOERR,IOLOG, ITER, 

X  ITER  1 ,  ITER2 ,  LENMA ,  LENHI ,  LENMY ,  LOOK ,  (•»,  MA  PA  (  20  ) , 

X  MAPI(20)  ,MAXA,MAXM,MAXN  ,N  ,NTYP£2,P,PRl4NT,TERMlN , 

X  UNBDDQiBASlS( 101 ) ,CAND(5) ,CANDI( 101 ,5) ,CANDL(5) , 

X  COLIC  101 ) ,ROWTYP( 101 ) ,STAT0S(331 ) 

C 

real  value 

INTEGER  10IN , NCOLSA , TT , SS 

COMMON  STIMLl , KW , WEIGHT , RhALF , Z2 , HEALER , OMEGA , B , EASCB , BASLE 
- ,bASUb .BETAR ,CAKDA,CANi/CJ .COLA , MEMORY ,UZERO,XbZERO , YQ.MAPA , 

-hAP I , BASIS , LAND , CANDI , CaNDL , COLI , ROWTI P , STATUS , KWDEF 
C 


n  o 


C  WITH  ChAh.uc.ii  lt<  Thh  EELGW(XMP)  PAhAi-.bTERS  ThE  SIZE  OF  Thb  PROBLEM 

C  lb  LlhlTcD  Of.LY  Eli  The  AVAlLAblLITl  Or  COMPUTER  CORE. 

C  FOh  PitGLLbr.S  LARGER  Tt;b  30  KeoThlCTbL  ATTHlBuTbSC  Vj  UNRESTRICTED 
C  ATTHIEUTtS)  Ai.D  100  PAIRWISE  COMPAKlStONS  SET  Thb  BELOW  PARAMETERS 
C  TO  THE  FOLLOW  l  t.G  VkLcES. 

C  f-iAXA=fiAXi.TbiiiAXFPC-riv.nXh 

C  KAAM=hAXFPC+1 

C  F:AAN=MAXW7S+3fcMAXFPC-*-l 

C  COLMAX-hAXH 

C  LuuK=iiAXi. 

C  LEEMY>=1 1  *£-iAXF.+2*MAXh+2i't-1AXA+o+3i'  (  X  ''DbNSE,ihAXK**2  ) 

C  DbhSb=MAXA/  ( MAXt-i^I-lAkh  ) 

C 

c 

C  IMTIALIZL  VARIALES  FOR  XMP 

C 

IOlH=5 

I0EHH=S 

IOLOG=9 

KAXA=3331 

HAXll=101 

HAXh=331 

C0LHAX= 101 

P=3 

LOOb=331 
FACTuh= 1 0 
LEi.hiY=20?CU 
bhDTYPsI 
PMNT=0 
C 

C  ThE  OUTPUT  WRITTEN  FOR  Thb  PRINT  STATEMENTS  WITHIN  XMP  WILL  EE 
C  PLACED  It;  A  FILE  CALLED  TAPE9,  ThIS  FILE  WILL  BE  A  LOCAL  FILE 

C  WKEii  ThIS  PROGRAM  IS  RUN  FhOM  A  TERMINAL.  TO  READ  Tht  CONTENTS 

C  OF  THIS  FILE  YOU  DO  THE  FOLLOWING  AFTER  ThE  PROGRAM  hAS  FINISHED 

C  RUNNING. 

C  REWIhD,TAP>  9 

C  COPY.TAPE9 

C  ThIS  FILE  IS  LOST  IF  NOT  MSA VE"-D  AFTER  ThE  RUN. 

C 

C 

C  INITALIZE  ALL  PARAMETERS  OF  ThE  MODEL 
HAXATT=30 
MAXSTI=100 
MAXFPCs 1 CO 
HAXWTS=30 

C  ZERO  OUT  THE  KWDEF  AND  THE  WEIGHT  ARRAYS 

DO  7  KARRAY=1 .MAXATT  * 

KWDEF(KARRAY)sO 
WEIGHT { KARR AY ) =0 . 

7  CONTINUE 


VARIABLE  LIST 


■ih%i  t*MVr  YtfYi 


w  ititii*Witirltfcit>il>kaK»ir*>ti>»fctt<>aii>*i<>itfcaliaaw<r»i>»aw4»»*>tii»>wiii»i>i»i>it*lrit  akalritkit 

a 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


**  STlhLI(t , J)=S  ATTttlbUTE  J  WITH /RESPECT  TO  STIMULUS  I  HAS  A  SCORE 

*  Or  5. 

k  OHtUA(I,1  )=L  STIMULUS  L  IS  PRfc.rfc.RRED  TO  STIMULUS  K  IK  A  FORCED 

6  OhLGA(l,2)=M  PAIRED  COMPARISON. 

* 

*  w£iGhT(i)=T  Value  of  the  i'th  weight , gotten  FROM  XmP 


*  RhALr'(  I ,  J)=A 

* 

* 

*  Z2(I)=R 

» 

*  KWDEF(I)=K 

ir 

*  KATTA 

*  KSTkLI 

*  KFCONP 

*  NUhlnTS 

it 

* 

*  KWFLAG 

* 

*  Kiri  (  I ) 

it 

it 

*  KERROR 

*  KAXWTS 


COEFFICENT  Or  Tht  J'Th  V.SIGhT  V;ITh/RESPECT  TO 
THE  I  FORCED  PAIRED  COMPARISON.  I.E.  ThIS  IS  Tht 
COEFFICEtiT  OF  THE  J'Th  VARIABLE  AND  I’Th  CONSTRAINT 
VALUE  OR  AMOUNT  TnE  I'Th  FORCED  PAIRED  COMPARISON 
IS  NOT  SATISFIED  BX. 

K=1  , MEANS  THE  I'TH  WEIGHT  IS  UNDEFINED  K=0, MEANS  THE 
I'TH  WEIGHT  IS  DEFINED 
If  OF  ATTRIBUTES 
i>  OF  STIMULI 

¥  OF  FORCED  PAIRED  COMPARISONS 
V  OF  inhlGHTS  USED  IN  THE  L.P.(I.E.  WREN  A  WEIGHT 
IS  UNRESTRICTED  IT  WILL  TAKE  UP  TWO  OF  THE  NUMUTS  IN 
THE  L.P.) 

-1  ALL  WTS<=0  ,  2  ALL  WTS  UNRESTRICTED,  1  ALL  WTS>=0, 

0  ViTs  ARE  MIXED. 

THE  I'TH  WEIGhT  Is  OF  THE  A  PARTICULAR  DIRECTION  , 

WhICh  IS  OF  TnE  SAME  CODING  AS  FOR  KWFLAG (EXCEPT  FOR 
ThE  ZERO  VALUE) 

A  FLAG  USED  TO  CHECK  FOR  ERRORS  INT  ThE  INPUT  DATA  * 
MAX  NUMBER  Of  WEIGHTS  THAT  CAN  BE  USED  IK  SOLVING  * 


* 

*•  HAXSTI 

*  MAXFPC 

*  MAXATT 
**  hEADER 


THE  C.C.  MODEL  * 
MAX  NUMBER  OF  STIMULI  ThAT  CAN  bE  SPECIFIED  * 
MAX  NUMBER  OF  FORCED  PAIRED  COMPARISONS  * 
MAX  NUMBER  OF  ATTRIBUTES  ThAT  CAN  BE  SPECIFIED  * 
AN  ARRAI  TO  HOLD  BO  CHARACTERS  OF  PnOBLEH  DISCRETION  * 


¥¥¥¥¥¥¥¥**¥¥¥¥¥¥¥¥¥¥¥  wwkkit¥¥mii  ¥»¥¥¥¥¥¥ 


WRITE  OUT  THE  HEADING  FOR  ThIS  PROGRAM 


WftlTE(b,55) 
55  FORMAT(IHI) 
WRITE(6,215) 
215  FORM AT ("  ++> 


WRITE(6,130) 

130  FORMAT  ( 1h  ,"C  0  M  ?  0  S  I  T  E’,,10X,"  CRITERION'', 
-10X,"*  KOCH  M") 
l«RITE(b,12B6) 


m 


\ 


120 b  FORMAT ( 1  h  ,//// ,21X,"kftlTT Eli  bl:  h.  X.  RuNYAN",//, 

-2U, "FACULTY  SPONSOR:  L.  F.  MCGINNIS’* ,// , 

-21X,"G£GnuIA  INSTITUTE  OF  TECHNOLOGY",//, 

-2 IX, "VERSION  1 " , // ,2 IX , "NOV  tMB&K  1b,  I960  ",///) 

HRITE(b,213) 

C  R»:AL  IN  THE  HEADER  UR  PROBLfcK  STATEMENT  OF  Th£  PROGRAM 
HEAD(3, 1740)  (HtAUER(lhEAD) ,IhEAD=1 , 10) 
mO  FORKAT(IOAo) 

READ(5, 10)  KATTA , KsTnLI , kFCOkP 
10  FOhhAT(3Ib) 

C  CHECK  FOR  ZERO  VALUES  OF  KATTA ,  KSTrtLl ,  KfCOKP 
IF(KATTA.Le.O)  GO  TU  44 
IF (KSThLI .LE .0 )  00  TO  44 
IF(KPCCMP.Lt.O)  GO  TO  44 
00  Tu  77 
44  inRIT£(6,6o) 

06  FOflhATdhC,"  ChECK  THE  VARIABLES ,  KATTii A ,  KSThLI , KFCOKP ,  FOH  VALUES 
-LESS  THAN  OR  EQUAL  TO  0.  -ERROR  S") 

KERR0h=1 
GO  TO  9999 
77  CONTINUE 
NUMWTS=KATTA 

READ( 5 ,13)  XwFLAG ,  (  KV. ( IsT  )  ,  IST=  1 ,  KATTA  ) 

15  F0RI-iAT(32I2) 

’  IF ( KLFLaG . £Q .0 )  GO  TO  1bb2  .  ^ -  % 

C  INITIALIZE  Th£  VALUES  OF  tACH  KV.(I)  IF  Kv.FLAG  IS  NOT  EQUAL  TO  0. 

DO  Ibbl  IST=1,  KATTA 
KV<  ( IST  )  skKFLAG 
16b1  CONTINUE 
1682  CONTINUE 
C 

C  ERROR  MESSAGES 

C  *  » 

C  *  ERROR  1  MORE  THAN  KAXATT  ATTRIBUTES  hAVE  BEEN  SPECIFIED  *> 

C  *  ERROR  2  MORE  THAN  hAXSTl  STIMULI  hAVE  BEEN  SPECIFIED  - 

C  *-  ERROR  3  MORE  THAN  MAXFPC  FORCE  PAIRED  COMPARISONS  hAVE  BEEN  SPECIFIED* 


C  *  ERROR  4  INVALID  VALUE  OF  KV.FLAG  hAS  BEEN  SPECIFIED  * 
C  *  ERROR  5  NUMV.TS  IS  CHEATER  ThAh  MAXWTS  * 
C  *  ERROR  6  A  SPECIFIED  HEIGHT  ,I.B.  KVi(I),  hAS  AN  INVALID  VALUE  • 
C  *.  ERROR  o  A  STIMULUS  IS  OUT  OF  RANGE  FOR  A  FORCED  COMPARISON  * 
C  *  ERROR  9  EITHER  KATTA, STIhLI, OR  KFCOMP  IS  <=0  * 


C 

C  BEGIN  CHECKING  FOR  OThER  ERRORS 
C 

KERRORsO 

IF (KATTA . LE .HAXATT )  GO  TO  20 
kRITE(b,30)  KATTA, KAXATT 

30  FORMAT ( IhO ,  16 , "  ATTRUIBUTES  HAVE  BEEN  SPECIFIED,  THIS  IS  IN  EXCESS 
-  OF  THE", 15,"  ATTRIBUTES  ThAT  CAM  BE  SPECIFIED.  "ERROR  1") 
KtRRORsI 

20  IF ( KSTHLI . LE . MAXST I )  GO  TO  33 


o  o 


WRITb(fa,45)  KSTMLI , NAXST l 

45  FORMAT ( 1  hO ,  lb , "  STIMULI  HAVE  fcbbt.  SPECIFIED,  THIS  IS  Hi  EXCESS  OF 
-The.",  16,"  STIHLT  TnAT  CAN  fab  SPECIFIED.  ‘•ERROR  2") 

KEKkORsI 

35  IF ( KFCOMP .  Lb  .ftAXr PC )  GO  Tu  50 
WnlTE(b,faO)  KFCO.'.P , rtAXr PC 

bO  F GhMAT (IhO.Io,"  FORCED  PAIRED  COMPARISONS  hAVb  faEEN  SPECIFIED,  THI 
-S  IS  IK  EXCESS  OF  Tnb  ALLOWAbLt.,,,Io,,,.  ‘'EhltOK  3") 

KLKHOR=1 

50  Ir (KWFLAG.EQ.-1 )  GO  TO  52 
IF  (  KWFLAU .  £<J .  0 )  GO  TO  52 
IF(KWFLAG.Ey.l)  GO  TO  52 
IF(KWFLAG.EQ.2)  GO  TO  52 
WhITE(6,110)  KWFLAG 

110  FORMAT (1hO,M  AN  INVALID  VALUE  OF  KLFLAG=,,,IV  HAS  faEEN  SPECIFIED. 

-  “ERROR  4") 

KEKR0ii=1 

52  CONTINUE 

IF(KERS0R.£Q.1)  GO  TO  S939 
IFUWFLAG.NE.O)  GO  TO  12b7 
DO  2231  OUKPs 1.KATTA 
KLZ=kW( JUMP) 

IF(KLZ.EQ.-I)  GC  TO  2231 
If(KLZ.EQ.2)  GO  TO  22 31 
IF(KLZ.EQ.I)  GO  TG  2231 
WHlTE(b,5431 )  JUMP ,KLZ 

5431  F0KAAT(1h  ,"AN  INVALID  VALUE  FOR  A  SPfcCIFIED  " 

-."WKIGhT  DIRECTION  HAS  faEEN  EhCOUnTEAEC",/ 

CttfcCK  THE" ,14,"  ENTER!  FOR  A  VALUE  OF", 13,"  “EHROft  6") 
XEhhORsI 
2231  CONTINUE 
1267  CONTINUE 

C  ChECK  TO  SEE  IF  ThE  LIMIT  HAXWTS  hAS  fcEEh  EXCEEDED 
IF(IABS(KWFLAG)-1)  fa 5,120,75 
C  NhEN  ALL  LEIGhTS  ARE  UNitbSTKlCTED 
75  L22shAXUTS/2 

IF(KATTA.LE.L22)  GO  TO  120 
GO  TO  65 
65  KKTEKPsO 

DO  90  ITT=1,NUMWTS 
KKTfcliP=KKTEMP+IAES(KW(  ITT) ) 

90  CONTINUE 

IF ( KKTfahP . LE . f-iAXWTS )  GO  TO  120 
C  THE  LIMIT  HAXWTS  hAS  BEEN  EXCEEDED 
65  V.R1TE(6,105) 

105  FURMATdKO,"  A  MAXIMUM  OF  30  WEIGHTS  CAN  BE  SPECIFIED  WHERE  EVERY 
-UNRESTRICTED  WEIGHT  COUNTS  AS  2."/"  THIS  30  WEIGHT  LIMIT  hAS  faEEN 

-  EXCEEDED.  ‘'ERROR  5") 

KERRGKsI 

120  IF(KERROR.LQ.I)  GO  TO  9999 

READ  IN  THE  SET  OF  FOhCED-PAlKbD  COMPARISONS (OMEGA)  CHECKING  FOR 
EKkOHS  AS  EACH  MEMBER  OF  THE  SET  IS  READ  IN. 

DO  140  IPKPsI ,KFCOmP 


o  n  o  o  on  n  n  r>  n  n  non 


Hc.fi.ul 5, 13b)  (Onega(IF.\P,1),I=1  ,2) 

1 3i>  F0rti-.AT(2I4) 

TT=0t‘iEGA(lPA.P,1) 
lf(TT.LT.I)  GO  TO  14b 
IF(TT.UT.KSTiiLI)  GO  To  145 
SS=GMEUA(  iFKF ,2 ) 

IF(SS.LT.I)  GO  TO  145 
IFISS.GT.kSTmH)  GO  TO  14b 
GO  TO  140 

K5  WrITe(u,  150)  1PKP ,OKbGA( IPkP ,  1 ) ,UH£UA( 1PAP ,2 ) 

150  FORMAT ( 1  h  1 , "  Thu.”,  14,"  FORCED  fAIF.tD  COMPARISON  CONTAINS  A  STlhULU 
-S  THAT  IS  OUT  Of  RANGE  CHECK  (“13,”  ,13,") .  *£RR0a  6") 

KEfinOH=1 

IF ( KbhftOR .  c.Q .  1 )  GO  TO  9999 
140  CONTINUE 

hEAD  IN  ALL  THE  VALUbS  FOH  ThE  ATTRIBUTES  ONE  STIMULUS  AT  A  TIKE 
DO  155  KIT=1 ,KSThLI 

READ(5,1bO)  (STIMLK  KIT , JEST ) ,JEST=1 ,KATTA) 
loO  FORMAT (bF 10.0) 

155  CONTINUE 
M=KFCGmP+1 

DETERMINE  THE  COaFFICENTS  OF  THE  WEIGhTS  IN  THE  CONSTRAINT  MATRIX 
El  SUBTRACTING  THE  SCORE  OF  A  CERTAIN  ATTRIBUTE  W/R  TO  A  STIMULUS  I 
FOH  THE  SCORE  OF  ThE  SAME  ATTrtlbUTb  V./R  TO  A  STIMULUS  J  WHEN  J  IS 
PHEFtrtED  TO  I  IN  A  FORCED  FAIRED  COMPARISON 
IF(IAbS(KWFLAG)-1)  1b5,170,1?5 
ALL  WEIGHTS  ARE  EITHER  NEGATIVE  OR  POSITIVE 
170  DO  100  IAbCs 1 , KFCOMP 
XFIRSTsOMEGA ( lAbC , 1 ) 

KSEChD=OMEG A ( IAbC , 2 ) 

DETERMINE  ThE  I'Ttt  COEFFICENT  ROW  OF  THE  CONSTRAINT  MATRIX  WHERE  ROW  I 
CORRESPONDS  TO  ThE  I'Th  FORCED  FAIRED  COMPARISON 
DO  100  JDEF=1 jNUMWTS 

RHALF(IAbC, JDEF)=XWFLAO"(STIMLI(KFIhST,  JDEF )-STIMLI(KSECND, JDEF)) 

KEEP  A  RUNNING  TOTAL  TO  DETERMINE  ThE  M'TH  ROW  OF  THE  CONSTRAINT  MATRIX 
RhALF ( M, JDEF )=RHALF ( H , JDEF )+RhALF ( IAbC , JDEF ) 

100  CONTINUE 
GO  TO  210 

ALL  WEIGHTS  ARE  UNRESTRICTED 
175  NUMWTS=2*NUKWTS 

DO  200  lABCs 1 , KFCOMP 
KFIRST=OMEGA ( IABC , 1 ) 

KSECKDsOMEGA ( IAbC , 2 ) 

DETERMINE  ThE  COEFFICENT  OF  ThE  CONSTRAINT  MATRIX  WHEN  ALL  ThE  WEIGHTS  ARE 
UNRESTRICTED 
DO  200  JDEFsl  ,NUMViTS,2 
C  COEFFICENT  OF  V.  EIGHT  (JDEF)+  AND  CONSTRAINT  IAbC 
KSCTs(JD£F+1)/2 

RhALF( IAEC , JDEF ) sSTlMLl ( KFIRST , KSCT ) -ST IMLI ( KSECUD , KSCT ) 

C  COEFFICENT  OF  WEIGHT (JDEF)-  AND  CONSTRAINT  IABC 
JDEF 2s JDEFsl 


oo  n  n  n 


RK  ALF ( IABC , JU&r 2 ) = -KhALF < IABC , JDEF ) 

C  A  hUI.filivC  TOTAL  Of  Tr.E  h'Th  EGG  OF  ThE  CONSTRAINT  MATRIX 
RHALF<(->,  JDEF)  =  RhALr  <H,  JDEF )+MhALr  (IABC  ,  JLtF) 

KhALF (M , J DEF2 ) =-KhALF< M , JDEF ) 

200  CONTINUE 
GO  TO  210 

C  THE  WhIGhTS  ARE  MIXED 
1b5  00  300  IABC=1,KFC0M? 

KFIRST=OhEGA( IAbC,  1 ) 

KSECNO=G.\iGA<  IAEC ,  2  ) 

KTfcMP=0 

C  KTEMP  IS  A  RUNNING  TOTAL  Ct  ThE  NUMBER  OF  COLUMNS  TAKEN  OP  b\  THE  WEIGHTS 
i>0  300  JDEF=1,NUi:WTS 
IF<IABS<KW< JDEF) ).CT.1)  GO  TO  1o5 
C  ThE  JDEF'TH  WEIGHT  IS  ElThth  NEGATIVE  OR  POSITIVE 
KT£KP=KT£fiP+1 

KhALF  <  IALC ,  KTEMP )  =KW  <  JDEF  )*(  STIMLI  <  KF  IhS? ,  JDEF  )  -STIHLI  <  KSBCND , 

- JDEF ) ) 

RHALF<M, KTEMP )=HhALF<h, KTEMP )+RHALF< IABC, KTEMP) 

GO  TO  300 

C  THE  JDEF'TH  WEIGHT  IS  UNRESTRICTED 
1b5  KTEHP=KTEMP+1 

RKALF< IAbC , KTEMP ) =S?lMLI(XFIhST , JOEf )-STIMLI<KSECND, JLEF ) 

RHALF  <  M ,  KTEMP  )  =  RHALF  (  H ,  KTcMP )  +HHALF  <  1 AEC ,  KTEMP  ) 

ISSTsKTEMP 

KTEMPsKTtMP+1 

KhALF < IAbC , KTtMP )= -KhALF < lAEC , ISST ) 

RH  ALF<  11 ,  KTEMP  ) = -KK  ALr  <  H ,  ISST  ) 

300  CONTINUE 

WUMWTSsKTEMP 
210  CONTINUE 

CHECK  TO  SEE  IF  TnE  h'Trt  ROW  OF  RhALF  IS  ZERO, IS  SO  ThEN  ThIS  MEANS 

THAT  THE  h'Th  CONSTRAINT  CANNOT  bE  SATISFIED  THEREFORE  THE  PRObLtti 

IS  iNFEASIbLE 

DO  1234  IChCKsI .NOMWTS 

IF<HhALF<M, IChCK) .NE.O)  GO  TO  1235 

1234  CONTINUE 
WRITE<6, 1212) 

1212  FCRMAT<1h  ,//,"  ThE  LEFT  HAND  SIDE  OF  ThE  M'TH  CONSTRAINT  IS  " 
-"ZERO.",/,"  ThE  STRUCTURE  OF  THE  DATA  SET  CAUSES  ThE  PROBLEM" 

-"  TO  BE  INFEASIBLE") 

.  GO  TO  9599 

1235  CONTINUE 

CHECK  TO  SEE  IF  THE  M'Th  ROW  OF  RHALF  CAN  SATISFY  THE 
RESTRICTION  ON  ThE  WEIGHTS, IF  ONE  EXIST 
IF<kWFLAG.EQ.2)  GO  TO  1999  * 

DO  2111  ICKCK=1 ,KATTA 
1F(KW< IChCK). £0.2)  GO  TO  1999 
IF<RhALF<M, IChCK ).GT.O)  GO  TO  1999 
2111  CONTINUE 

C  ThE  PROBLEM  IS  INFEASIBLE, IT  CANNOT  SATISFY  THE  H'Th  CONSTRAINT 
V.'RlTE<6,fc6o6) 

o Ciiti  FORMAT <  1h  ,//,"  BECAUSE  OF  THIS  PARTICULAR  SPECIFICATION" 


noon  on  oo  on  oooooo  n  o  oe 


-•»  UK  WEIGHT  RESTRICTIONS" ,  / , "  THE  PROBLEM  IS  INFeASIBlE" 

LETTIf.C  ALL  WEIChTS  Lb  UhhtSTHlCTtD  WILL  WORK") 

GO  TO  yyyy 

lyyy  continue 

SEND  ALL  ThE  DEFINED  COLUMNS  OF  THE  FIRST  r.UiiWTS  COLUMNS  OF  fthALF  TO  X«P 


CALL  AMAPS(  Hm/TYP ,  IOERh ,  LEt.MA ,  LEriHI ,  LENnl  ,NAPA  ,KAPI ,  MAX  A  ,MAAM 
-,IiAXN,MEMOhY) 

Tht  FIhST  K-1  CONSTRAINTS  Ahb  >=0  WITH  A  HhS  VALUE  OF  0. 

DO  MOO  11=1  ,kbCCi\r‘ 

KOWTYPlII}=-1 

b(ll)=0. 

MOO  CONTINUE 

ThE  M’TH  CONSTRAINT  IS  =1 

ROWTYP(M)=0 

B(M)=1. 

N=0 


DETERMINE  ThE  NUMBER  Cc  NONZERO  ELEMENTS  ( COLL  Eli )  IN  ThE  FIRST  NUMWTS 
COLUMNS  OF  ThE  CONSTRAINT  MATRIX*  ALSO  DETERMINE  WHICH  HOWS(I.£.COLI)  HAVE 
NONZERO  ELEMENTS  IN  THEM  FOR  A  CERTAIN  COLUMN  AND  WHAT  THAT  NONZERO 
ELEMENT  IS(I.E.  COLA) 

CJ=0. 

CJ  IS  ThE  OBJECTIVE  FUNCTION  COEFFICEf.T 
DO  600  JUST=1 .NUmVTS 
CULLEN =0 
DO  500  Lt'iNsI  ,M 

IF( RhALF( LmN , JUST ) . EQ .0 )  GO  TO  500 

COLLEMsCCLLfcN+1 

COLA ( COLLtN )= R HALF ( LMN , JUST ) 

CCL1  (  COLLEN  )  =  LJJ* 

500  CONTINUE 

CHECK  FOR  TO  SEE  It  THIS  COLUMN  hAS  ALL  ZEROS  IN  IT 

ALSO  FLAG  ThE  UNDEFINED  WEIGHTS 

IF(COLLEN .EQ.O)  KWDEF< JUST)=1 

IFCCOLLEN.BW.O)  CO  TO  600 

PASS  THIS  COLUMN  TO  AMP 


CALL  XADDAJ ( C J , COLA , COLI , CULLEN , COLMAX , IOE  Rfi , J , LENNA , LENMY , MAPA , 
-MEMORY, K) 


600  CONTINUE 

CHECK  TO  SEE  IF  THERE  EXIST  ANY  WEIGHTS  THAT  ARE  DEFINED 
IF(N'.NE.O)  GO  TO  602 
WRITE (6, 601) 

601  FORMAT ( 1h  ,  "ALL  V.EIGhTS  ARE  UNDEFINED  V.hICH  MAKES  IT  IMPOSSIBLE* 

TO  SOLVE  THIS  PARTICULAR  PROBLEM  WITH  THIS  DATA  SET*) 

GO  TO  9555 
002  CONTINUE 

SEND  THE  REMAINING  COLUhNS(I.E.  THE  Z'S)  TO  XMP 

NOTE:  EACH  Z  COLUMN  HAS  ONE  NON-ZERO  ELEMENT  AND  AN  OBJECTIVE  FUNCTION 
COEFFICENT  OF  1.EUT  SINCE  XMP  ONLY  MAXIMIZES  WE  MUST  USE  -1  AS  THE 
OBJECTIVE  FUNCTION  COEFFICENT 


ooon  o  o  coono 


C0LLEN=1 

CJ=-1. 

C0LA(1)=1. 

BO  700  J2=1,KFC0riP 
CGL1(1)=J2 


CALL  XADDAJ  ( C J , COLA ,  COLI , COLL Et, , COLMAX ,  IOEkR  , J ,  LENMA , LENHY , 
-MAPA.MEMOhY.N) 

C  - 

700  CONTINUE 


DO  710  JXJ=1,N 
STATUS(JXJ)=0 
710  CONTINUE 


CALL  XSLACK ( B , E  A3CB , BASIS , BASLE , B ASUb , BNDT YP , BODED , COLA , COLI , 

X  COLHAX , IOERH ,  LENMA ,  LEI. HI ,  LEfcfcl ,  M ,  MAP  A  ,MAPI ,  HAXM , MAXK ,  MEMORY ,  H , 
X  KOViTYP ,  STATUS ,  UZEhO ,  XB2ERG ,  Z  ) 


SOLVE  THE  L.P.  PROBLEM 


CALL  XPKlML(b,£ASCf>, BASIS, BASLb.BASUB ,BNDTYP,BOUND,CAND,CANDA, 
X  CANOCJ.CAf.DI  ,CANLL ,  COLA ,  COLI  t  CGLMAX ,  FACTOh ,  IOEhR  ,  IOLOG ,  ITER1 , 
X  ITER2 ,  LENMA ,  LEhKI ,  LEf.Mi ,  LOCK ,  h.MAPA ,  MAPI  ,MAXm ,  MAXK ,  MEMORY  ,h , 
X  NT  YPE2 ,  P ,  PhlKT ,  STATUS ,  TEhMlN ,  ONBBUi ,  UZERO  .XbZERO ,  IQ ,  2  ) 


bRITE(6 ,7*42 )  (HEADER <  IhEAD ) , IHEAD=  1,10) 

7H2  FUhMAT(1h1 ,// ,  IX, IGAo ,////) 

WRITE  OUT  THE  TERMINATION  CODE  OF  TtiE  L.P.  PRObLEfi  TO  TAPE  S 
WHITE (9, 720)  TERilIN 

720  FORMAT ( 1  h  ,  "THE  TERMINATION  CODE  IS", 15) 
lF(TEmiIN.EQ.1 )  GO  TO  7*i0 
IF(T£RMIN .HE .4 )  GO  TO  S995 
WRITE(6,135&) 

1358  F0RMAT(1H  ,"***WAftNING*fc*",10X,"fc**WARNING*fcfc",10X,"***WARNING«»«" 
1 0X,  "^•’WARNING*"*" ,/, "  THE  OPTIMAL  SOLUTION  DOES  NOT  SATISFY", 
THE  ACCUHACY  CHECK  FROM  XMP , TERMINATION  CODE  4.") 


TEfihlN  IS  THE  TERMINATION  CODE  AND  IS  EXPLAINED  IN  THE  DOCUMENTATION 
OF  XMP 


C  DETERMINE  ThE  VALUES  OF  THE  HEIGHTS  FROM  ThE  XMP  VAKIAbLES 
740  IF(IAbS(KNFLAG)-1)  725,730,735 

C  . . . 

C  KWT  KEEPS  a  RUNNING  TOTAL  of  what  xnp  variable  YOU  are  TALKING  AEOUT 
C  WITH/RESPECT  TO  THE  ARRAY  STATUS. 

730  KfeTsO 

C  ALL  WEIGHT  ARE  EITHER  NEGATIVE  OR  POSITIVE 
DO  BOO  1JJ=1 .UUMWTS 

C  CHECK  TO  SEE  IF  ThE  IJJ’Th  hEIGhT  IS  UNDEFINED 


n  a 


IF<KWDfcF(IJJ).fcQ.1)  UO  TO  600 
Kl;T=Kl.'T+1 

C  CHECK  TO  SEE  It*  The.  IJJ  VaRIALLE  IS  IK  Thu  BASIS 
IF ( STATU S Ul.T ) . fcIC, . 0  )  uO  TO  7^5 
lr.  ElUhT  ( 1 J  J  ) sFLUAT  (  KV.FLXG  )  *XbZbi(U  (  STVi'uS(  KViT ) ) 

UO  TO  600 
74 5  WElUHT(lJJ)=0. 

O00  CONTINUE 
GO  TO  GOl 

c  . 

C  ALL  WEIGHTS  A. HE  DNhESTKlCTbD 
735  KOuftT= 1 
KVjT=0 

DO  1000  IJJ=1,N0hWTS,2 

C  CHECK  TO  SE£  If  THE  hGuNT'Tf.  WElUhT  IS  LNDEFlhtb 
It  (K'»DEF(IJJ  )  .EQ.  1 )  KWDEF(KOUNT)=  1 
IF(KNBEF(IJJ).£Q.1)  GO  TO  757 
KWT=KWT+1 

C  ThE  KGUhT '  TH  V, EIGHT  IS  DEFINED 
KWDEF(kCONT)=G 

C  DETERMINE  TfcE  XGUNT'Tn  WEIGHT  IF  IT  IS  DEFlKe.0 
LlKskWT 

765  IF( STATUS ( KV.T ) . EG . 0  )  GO  TO  7*47 
C  THE  Hub  NT' Th  V,t.IGHT  IS  NONZERO 
W  EIU  HT  (  KWT  )  s  XE2  ££0  (  ST  aTus  ( «i.WT  )  ) 

GO  TO  760 
7**7  KElGhT(EWT)=0. 

750  CONTINUE 

CHECK  TO  SEE  IF  WE  HAVE  ALREADY  CHECKED  £0?h  ThE  POSITIVE  AND  NEGATIVE 
PAHTS  OF  THE  KGUNT'Th  WEIGHT 
IF( LIK.NE . KV.T )  GO  TO  755 
KWT=kV»T+1 

C  WHEN  LIk=KV.T  W£  ARE  CHECKING  TO  SEE  IF  ThE  NONNEGATIVE  PART  OF  THE 
C  KODriT'TN  WEIGHT  IS  POSITIVE, SIMILARLY  IF  KWT=L1K+1  WE  ARE  CHECKING  TO 

C  SEE  IF  ThE  NON POSITIVE  FART  OF  THE  KOUNT'Th  V.EIGhT  IS  NEGATIVE 

CO  TO  765 

755  WEIGhT ( KObKT )= WEIGHT ( LIK ) -WEIGHT ( KWT ) 

757  K0UN?=K0DNT+1 
1000  CONTINUE 
CO  TO  oOl 

C  . 

C  ThE  WElGhTS  ARE  MIXED 

725  KT£iiP=1  .  .  . 

KFLAGsO 

C  KTEhP  KEEPS  A  RUNNING  TOTAL  Ori  ThE  NUMBER  OF  WEIGHTS, ACTUAL  WEIGHTS  AND 

C  FOR  EXAMPLE  IF  Th£R£  ARE  10  WElGhTS>=0  Ar.D  5  DNREStlCTED  WEIGHTS  ThEN 

C  ThE  FINAL  VALUE  OF  KTEMP  WILL  BE  15. 

KWT=0 

DO  1100  MN0=1,U0hWTS 
1F( KFLAG . NE . 1 )  UO  TO  1730 
KFLAGsO 
GO  TO  1100 

1730  IF(IABS(KV.UTtMP)).EC,.1)  00  TO  1745 


o  o 


C  Thb  KTEMP'TK  WhIGhT  Is  UKNLSTmICTEO 
KFLAG=  1 

C  CHECK.  TO  SEE  IF  ThE  KUuM'Tn  HtlGtiT  15  UNDEFINED 
lF(hWDEF(MNO).EG.1 >  KWDbFUuDJ.T)=1 
It(KW0fiF(M*\0).bQ.1)  uu  Tu  1770 
KwT=KWT+1 

C  Tub  kOUNT'Tii  WEIGHT  IS  DEFINED 
KUCEF  (  KOUNT  )  =0 

lFtSTATuSU»«T).EC.G)  GO  TO  1765 
WElGHT(K?£Kp)=Xb2ahU(sTATuS(XWT)) 

KY.T=KV,T+1 
GO  TO  1770 
1755  KWT=KWT+1 

IF ( STATUS ( NUT ) .  HE . 0 )  GO  TO  1766 
WEIGHT ( KTEMP ) =0 . 

GO  TO  1770 

1766  WEIGHT (KTEinP  )=-XEZEhG  (STATUS(  KV.T )  ) 

GO  TO  1770 
17^5  CONTINUE 

THE  XTEMP’Trt  WEIGHT  IS  EITHER  -  OH  +  ,BUT  IS  HOT  UNRESTRICTED 
CHECK  TO  SEE  IF  THE  XGGM'Th  WEIGHT  IS  UNDEFINED 
IF(KUDSF(mNO)  .EQ.1 )  XUDEr  U0UNT)=1 
IF(KWDEF(hN0).bCi.1)  GO  TO  1770 
KWT=KWT+1 

IF  (  STATUS(  XWT ) . HE . 0 )  GO  TO  1775 
VibIUhT(KT£tfP  ) =0  . 

GO  TO  1770 

1775  KEUihT(KTLS? )  =KVj  ( KTEEP ) -XfaZbhu ( STATUS ( XWT ) ) 

1770  KT£t-iP=KTEhF+1 
1100  CONTINUE 

C  . . 

601  CONTINUE 

C  OETbKrilNa  ThE  VALUE  CF  THE  Z'S  FOh  ALL  KViFLAGS 
KSCRFTsNUhwTS 
DO  22  IJK=1,KATTA 

IF  (  KWDbF  ( I JK  ) .  EG .  1 )  KSCHPTsKSCliPT-  ( IABS(  KW ( IJK ) ) ) 

22  CONTINUE 

DO  500  IIJs1,KFC0:-P 
KSCRPT=KSCRPT+1 

IF(STATUS(KSChPT) .EG.O)  GO  TO  615 
Z2 ( XI J ) =XbZERO ( STATUS ( KSCHPT ) ) 

GO  TO  900 
615  Z2(IIJ)=0. 

900  CONTINUE 

C  SINCE  WE  USED  XMP  TO  MAXIMIZE  AND  IN  REALITI  WE  WANTED  TO  MINIMIZE  WE  HAD 

C  TO  MULTI PI  THE  OBJECTIVE  FUNCTION  COEFFICENTS  bl  -1  SO  OUR  Z(OPTIMAL) 

C  OBJECTIVE  FUNCTION  VALUE)  Is  hEALLX  =-Z 
Za-Z 

WP.ITE(6,1920)  Z 

1920  FORMAT { 1h  ,"ThE  OPTItiAL  OBJECTIVE  FUNCTION  VALUE  IS  Z=W,F10.5) 
t  WHITE  OUT  ThE  HEADING  FOH  ThE  OUTPUT 
WR1TE(6,1530) 

1530  FOhi-iATdh  ,//,10Xl,%EIGhT^^••,17X,,,**-VALUE^•«•,) 


khlTb(b,1<o5) 

1*36  funt.AT(  IK  ,9X," - ".ITa," - ") 

C  UOT  THE  V'ALUt  Of  Tht.  WhlGhTS 

C  IDtt’LU  IS  A  FLnU  WhlCr.  IS  hQUAL  To  uNE  IF  Thtt\E  AliE  A t«Y  ONDErlhbD  WtlGhTS 
lDct'LC=0 

DO  1940  IIJ=1,kATTA 

C  Chc.CE  TO  Stt  IF  Tht.  *\A?TA*Th  V.ElUhT  lb  UNbSFlfvtU 
lF(M.DFF(IIJ).tQ.O)  GO  To  193b 
WRIT£(6,  1*42)  11J 

1942  format O  sx ,  13 ,  i  ox ,  unde?  in  el  " ) 

IDtt'LG=1 
CO  TO  194c 

1936  WRITE (6,1*50)  IIJ ,WElGh?(lIJ) 

1950  FORMAT ( 1 2X , 13 , 1 ?X , F 1 0 . 5 ) 

1940  CONTINUE 

ir;hITE(6, 1*oC) 

1960  F0hNAT(1h1,////,l4k,M2,,,19Xt"*»»VALU&fc-") 

VtRITr.(6 , 19o1  ) 

1961  FORMAT (  1h  ,12X," - ".loX," - ") 

C  LRIT&  06?  Tht  VALUES  of  Tr.E  2'S 

DO  1200  IJKLM=1  ,ErCO.\P 
V.'RITE(6, 1165)  IJKLM,L2(  iJALM) 

1165  FORMAT ( 1 H  ,  1 1X,  13, 17a , MG .5) 

1200  CONTINUE 

C  CALCULATE  ThS  SCOAtS  ru.i  EACh  STIMULUS  AND  V.hlTE  THEM  OUT 
WHIT£(6,2J22) 

2022  FORMAT ( 1  h  1 ,////,SX, "STIMULUS" ,  m/'^SCOKb**") 

WRITt(o,2C23) 

2023  FORMAT ( 1  h  ,bX," - ",1?X," - ") 

DO  1202  KIT=1  ,KS?i.LI 

SCOkfc=0 . 

DO  1203  KRR=1,KATTA 

IF ( EWDEF (ERR).EQ.I)  GO  TO  1203 

SCOit  t = SCOA  E+N  EIGhT (  KRH )  ‘STIHLI  ( KIT , KHR ) 

1203  CONTINUE 
khITE(6,2229)  KIT, SCORE 

2229  FOHMATOh  ,11X,I3,1?X,f  10.5) 

1202  CONTINUE 

C  WRITE  OUT  A  MESSAGE  If  ANY  Or  THE  ViElUhTS  ARE  UNDEFINED, I. E.  IDEFLG= 1 
IF(IDEf LG.EQ.O)  GO  TO  1204 
WRIT£(b,17) 

17  FORMAT(1h  ,//,"  FOR  Tht  DATA  SET  USED  IN  THIS  PROBLEM  SOME  OF  THE" 

ABOVE  ViEIGhTS  COULD  NOT  BE  DETERMINED") 

1204  ViRIT£(6,1976) 

1976  FORMAT ( 1 h 1 , /// , "  PROGRAM  TERMINATES  NORMALLY") 

5599  CONTINUE 
STOP 
END 
*  END 


